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ABSTRACT 



This thesis deals with the enhancement of video images 
degraded by turbid water viewing conditions. An algorithm by Peli 
and Lim has been used with some success for enhancement, but it 
was found to accentuate noise. The thesis examines a combination 
of the Peli and Lim algorithm with three approaches to 
enhancement. 

First, a Short Space Spectral Subtraction algorithm which 
performs the restoration in the density domain, using an estimate 
for the power spectrum of the given data set. The degraded image 
is divided into many subimages and each subimage is restored 
separately and then combined. 

Next, an algorithm for Image Enhancement and Noise Filtering 
by Use of Local Statistics, which uses the assumption that the 
sample mean and variance of a pixel is equal to the local mean 
and variance of all pixels within a fixed range surrounding it. 

Finally, a median filter for noise reduction , where a given 
pixel of a degraded image is replaced by the median of the pixel 
values in a window surrounding it. 

Combination of the algorithms are applied to degraded images, 
and the results are compared and discussed, in each case. It was 
found that noise smoothing can be achieved with the spectral 
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subtraction algorithm, and that the local statistics technique 
yielded very good contrast enhancement. 
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I. INTRODUCTION 



In this thesis, some algorithms for the enhancement of video 
images degraded by turbid water viewing conditions are 
implemented. The images used to illustrate the processes, were 
recorded during the recovery procedure in a torpedo testing range. 

The recovery operation is carried on undersea with the aid of 
digging equipment and underwater video cameras equipped with 
strong artificial lights. The equipment is controlled and the 
operation is monitored on video monitors at the surface onboard 
the recovery vessel. The recovery equipment in its attempt to dig 
out the torpedo, stirs up sediment wich visually obscures the 
object of interest and impedes the operation. To date no special 
techniques have been used operationally to process or enhance the 
video image before display. 

When it is of interest to “improve” an image, there are two 
broad types of image manipulation processes that cover all 
operations performed to get such improvement: image restoration 
and image enhancement. 

The goal of image enhancement is to process a degraded image 
so that the result is more suitable than the original image for a 
specific application or aids the human analyst in the extraction 
and interpretation of pictorial information [Ref. l]. Image 
restoration, on the other hand, is intended to bring a degraded 
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image back to an ideal degradation-free image as closely as possible 
[Ref. 2]. This makes the image enhancement problem a 
subjective one, the improvement of the image appearance to the 
human viewer is highly dependent on the viewer himself and on 
the application. What is “good” for one person or application is not 
necessarily “good” for another one. Consequently, the suitability of 
the processed image for a specific application makes image 
enhancement techniques very much problem-oriented. 

The degradation process can result from causes such as, an 
imperfect photographic process, imperfect display devices, poor 
contrast due to environmental conditions, different forms of noise 
(channel, quantization, salt and pepper) and others [Ref. 3], In 
the same manner, techniques for enhancement are dependent of 
the original degradation process. Processes yielding satisfactory 
results for one type of degradation are not necessarily suitable for 
another. 

The interest of this thesis is on examining different techniques 
to enhance images that are degraded by environmental conditions 
that create a great lack of contrast in the recorded images. This 
work is a continuation of previous research begun on the subject, 
which included the implementation of an adaptive filtering [Ref. 4] 
algorithm for contrast enhancement. This algorithm was applied to 
a set of underwater images recorded in turbid water viewing 
conditions. The images are low contrast degraded by noisy 
background. As will be explained in the next chapter, the 
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algorithm yielded very good contrast enhancement but, seemed to 
accentuate the noise. The purpose of this work is to investigate 
some noise smoothing techniques to be applied to the 
contrast-enhanced images, an alternative contrast manipulation 
scheme, and to compare the results obtained to those found with 
the algorithm already implemented. Besides the images included in 
the experimental results of this thesis, many more images 
resulting from other combinations and variations of the processes 
and different settings of the parameters were obtained. We have 
included here only those thought to be most representative of the 
results obtained. 
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II. OVERVIEW OF PREVIOUS WORK 



A. DESCRIPTION OF THE ADAPTIVE FILTERING ALGORITHM 

An adaptive filtering algorithm for image enhancement by 
Peli and Lim [Ref. 4] has been used for contrast enhancement of 
images. The images are typifyed by a large dynamic range, and 
recorded over a medium, characterized by a much smaller 
dynamic range, which causes those regions of the images with 
very high or very low luminance to be poorly represented. This 
algorithm modifies the local luminance mean of an image and 
controls the local contrast as a function of the local luminance 
mean of the image. 




Figure 2.1. Adaptive Filtering for Image Enhancement 



In Figure 2.1 a block diagram of the Adaptive Filtering 
algorithm is shown In the figure, f(nj,n 2 ) denotes the 
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unprocessed digital image, and fL(n 1 ,n 2 ), which denotes the local 
luminance mean of f(n!,n 2 ), is obtained by low-pass filtering 
f(n 1 ,n 2 ). The sequence fn(ni,n 2 ) which denotes the local contrast 
is obtained by subtracting fi.(ni,n 2 ) from f(nx,n 2 ). The local 
contrast is modified by multiplying fn(ni,n 2 ) with k(f L ), a scalar 
which is a function of f L (ni,n 2 ). The modified contrast is denoted 
by f'H(ni,n 2 ). The specific functional form of k(fL) depends on the 
particular application under consideration; a value of k(f^ > 1 
represents local contrast increase while k(fiJ < 1 represents local 
contrast decrease. The local luminance mean is modified by a point 
non-linearity and the modified local luminance mean is denoted by 
f' L (n 1 ,n 2 ). The specific non-linearity chosen depends on the 

particular application under consideration, and in most application 
problems the non-linearity is chosen so that the overall dynamic 
range of the resulting image is approximately the same as the 
dynamic range of the recording medium. The modified local 
contrast and local luminance mean are then combined to obtain 
g(n 1( n 2 ), the processed image. 

B. DISCUSSION OF RESULTS 

In this section the application of the algorithm described in the 
previous section to the images shown in Figure 2.2 is presented. 
Figure 2.2 shows two low-contrast noise-degraded images of 512 X 
512 pixels with each pixel represented by 8 bits. Note that the 
small dynamic range of the luminance of both images make 
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(a) 




(b) 

Figure 2.2 Original Degraded Images. 
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difficult the appreciation of subtle details in them In this case, it 
is desired to increase the local contrast in the images to bring back 
the lost details. To achieve this, k(f L ) is chosen as shown in Figure 
2.3 (a) and (b) for the images of Figure 2.2(a) and (b), 

respectively. For the case of Figure 2.2(a), the function is chosen 
to increase the local contrast for the low luminance regions, 
maintaining it relatively constant for medium ranges of local 
luminance, and again enhance it, for the brighter regions of the 
degraded image. The results are shown in Figure 2.4(a), where a 
relative contrast increase can be observed. This is especially 
evident in the “fish” shape close to the lower right corner of the 
image. Also note the increase in background noise resulting in the 
textured background observed in the contrast-enhanced image. 





(a) (b) 



Figure 2.3. K(f L ) Function for Processing Images in Figure 2.2 
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(a) 




(b) 

Figure 2.4. Images in Figure 2.2 Processed with Adaptive Filtering 
for Image Enhancement Algorithm. 
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In the case of the image of Figure 2.2(b), it is desired to 
retain the original local contrast in the darker regions, but 
enhance it in the regions corresponding to medium and high local 
luminance means. This explains the particular shape of the k(f L ) 
function chosen to process this image. In the results, shown in 
Figure 2. 4(b), some details difficult to be distinguished in the 
un-processed image, are now more evident. Among them are the 
particular shape of the object near the center of the image and 
the form of the ocean floor below the left-most object. Again, it 
can be observed that there is an increase in the background noise 
throughout the contrast-enhanced image. 
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Ill DESCRIPTION OF THE ALGORITHMS 



In this chapter, the theorethical foundations and procedures 
for implementation of the different algorithms used in the 
noise-filtering and contrast enhancement problem are presented. 
We start with a description of noise-smoothing techniques to be 
used as post-processing after the adaptive filtering algorithm. Then 
an alternative method for the contrast enhancement is presented. 
Finally, the use of median filtering for noise-smoothing is 
described 

A. SHORT SPACE SPECTRAL SUBTRACTION 

This algorithm, developed by J. S. Lim [Ref . 5] performs the 
restoration in the frequency domain, using a converging solution 
for the power spectrum of a given data set that does not depend 
on the original estimate. The procedure is based on two basic 
assumptions. First, each part of an image generally differs 
sufficiently from other parts so that the image cannot be modeled 
by a stationary random process. Second, the power spectrum of 
the restored image is estimated by the spectral subtraction of the 
additive noise spectra from the degraded image. 

1 . Spectral Subtraction 

Frequency domain techniques for image enhancement use 
iterative procedures for estimating the image power spectral 
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density, continuing until certain fidelity criteria are met The 
density so obtained, is dependent on both the initial estimate for 
starting the iterative procedure and the stopping criterion. Some 
methods converge more rapidly than others, and sometimes first 
and second order statistics of the degrading function must be 
known . 

Consider a model of image degradation as presented in 
Figure 3 . 1 . In this model f(n 1 ; n2) represents a digital image, 
b(nj,n2) represents a linear space-invariant point spread function, 
and d(n 1 ; n2) represents an additive noise component. Thus, 
g(ni,n2) represents a blurred version of f(n 1 ,n 2 ) and y(n 1 ; n2) is 
the degraded image. 




An image restoration system is shown in Figure 3 . 2 . It 
can be seen that the goal of a restoration system is to process a 
degraded image y(n 1 ,n 2 ) through a noise reduction system to 
estimate g(ni,n2) which in turn is processed by a deblurring 
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system to get f(n 1 ,n2), an estimate of the original blur-free noise- 
free image 




Figure 3 . 2 . Image Restoration System 



The purpose of this algorithm is to provide a model for 
the Noise Reduction System, as presented in the Figure. A general 
model for the restoration filter is given by 



H(« ,W 2 ) 



-{ 



P (CO ,(0 ) 
£_ 1 2 _ 



P (co ,c o ) + a • P (co ,co ) 



} 



( 3 . 1 ) 



1 2 



1 2 ' 



where P g (co^ , CO2) and Pdfai.^) represent the power spectra of 
g(ni,n2) and d(n 1 ,n2) respectively, a and P are constants. If P is 
unity, H (coj, C02) corresponds to the parametric Wiener filter, and if 
a also is unity, it reduces to the standard Wiener filter [Ref. 1 , 2 ]. 
Other values for a and P define different filtering techniques, such 

as power spectrum filtering and geometrical mean filtering. One 
commonly used approach in the implementation of H (00 1,002) is an 

iterative procedure which begins with an initial estimate of 
Pg (001,002) and then iteratively estimates g(ni,n2) and P g (ooi,oo2) 



18 



until a converging solution or a desirable performance is achieved 
This iterative procedure may be computationally undesirable. 

The Spatial Subtraction approach analytically obtains a 
converging solution, by estimating P g (001,002) from g(ni,n 2 ) as 
(i/k).|G(oo 1 ,oo 2 )I 2 and using a value for (3 of V2. The complete 
solution for this problem, is given by 



and zero otherwise 

Here G^,^) represents the discrete space Fourier Transform of 
g(ni,n 2 ), Y (ci)i,(i) 2 ) is the discrete space Fourier Transform of 
y(n!,n 2 ), represents the power spectrum of d(ni,n 2 ), 

the additive noise component of the degradation model, and k is a 
scaling factor that normalizes the power and energy spectral 
densities. The phase of g(ni,n 2 ) is estimated by the phase of 
y(ni,n 2 ) and the transform magnitude of g(ni,n 2 ) is estimated by 
a particular form of spectral subtraction. 




( 3 . 2 ) 



* G(w ,W) = * Y(U,W ) 



( 3 . 3 ) 



for 



|Y(w ,(0 2 )| 2 ;> a.k.P d (w .w) 
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2. Short Space Implementation 

As mentioned before, the basic assumption is that each 
part of an image f(ni,n2) generally differs sufficiently from other 

parts so that it cannot be modelled by a stationary random field. 
Thus a short space implementation seems appropiate. The 
technique consists of dividing the original image into many 
subimages and restoring each subimage separately. The procedure 
involves the application of a short space window function 
WjjOii,^) to overlapping portions of the degraded image yO^,^) 

so that: 

y(iVn 2 )«w M (n r n 2 ) = {gfn^n^+dfn^)} • ^(n^) (3.4) 
or, equivalently, 

= 9 i,;|K' n 2 > +d i,;|( n i' n 2 ) < 3 ' 5 > 

The noise reduction system is then applied to yij(ni,n2) 
to recover gijCni,^). The overall full-size restored image g(ni,n2) 
is given by 



2K 2L 

g(n r n 2 ) =22 # i( n i' n 2> 

i=o 3=0 



(3.6) 
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The window function is desired to be a smooth function to 
avoid some possible discontinuities that may appear at the 
subimage boundaries in the processed image and must satisfy: 

2K 2K 

2 2 w i i( n i' n 2 ) = 1 ( 3 - 7 ) 

i-0 j-0 

for all of interest 




The condition in Equation 3 . 6 guarantees that an image 
can be reconstructed from its subimages. To ensure the smoothness 
of the window, a 2-D separable triangular window similar to that 
shown in Figure 3.3 will be used. Each window overlaps its 
neighboring window by half the window length in each dimension. 
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B. CONTRAST ENHANCEMENT BY USE OF LOCAL STATISTICS 

This algorithm presents a computational technique for contrast 
enhancement on a two-dimensional image array based on the local 
mean and variance. The algorithm, as developed by J. Lee [Ref. 
6] uses the basic assumption that the sample mean and variance 
of a pixel is equal to the local mean and variance of all pixels 
within a fixed range surrounding it. 

Let Xij be the brightness of a pixel (i, j) in a two dimensional 

NxN image. The local mean and variance are calculated over a 
(2n +l)x(2m +l) window. The local mean is defined as 



The algorithm is designed such that a pixel Xjj will maintain 
its local mean, and yet permit its variance to be modified by a 
constant factor times its original variance. The procedure is defined 
by 




k-i-n 1-j-m 



(3.8) 



and the local variance is defined as 




(3.9) 




(3.10) 



22 



where k, the gain, is the ratio of new local standard deviation to 
the original standard deviation. This approach has the 
computational advantage that the local variance v 4 j is not 
required and only the local mean mjj needs to be computed. If 
we have k > 1, the image will be sharpened as if acted upon by a 
high-pass filter. If 0 < k < 1, the image will be smoothed as if 
passed through a low-pass filter. In the extreme case, one has k = 
0 and x'jj is equal to its local mean mjj . 




Figure 3.4. Contrast Enhancement by Use of Local Statistics 



Figure 3.4 shows a block diagram of the algorithm. The local 
mean is evaluated over a (2n+l) (2m+l) window and subtracted 
from the original brightness of the pixel, Xij. The difference is 
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then multiplied by the standard deviation ratio to enhance the 
local contrast. Finally, this result is added to the local mean 
yielding the reconstructed pixel value x'jj. 

C. NOISE FILTERING BY USE OF LOCAL STATISTICS 

This algorithm, also developed by J.Lee [Ref. 6] is actually a 
variation of the previous algorithm adapted for the noise filtering 
problem. In this algorithm, the a priori mean (or variance) of 
the estimated image is calculated as the difference between the 
mean (or variance) of the noise corrupted image and the mean 
(or variance) of the noise itself. Let Zjj be the degraded pixel x i} j. 

Then the degraded pixel can be modelled as the sum of the 
noise-free pixel and a white noise sequence. That is, 



where wjj is the white random sequence with Efwjj ]=0. Define 
the a priori mean and variance of xjj to be 




(3.12) 



and 



°i.r e i < 2 




(3.13) 



respectively. The estimated pixel value Xjj is computed by 
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( 3 - 14 ) 



; i.r x "i,j * k i.] (z i.r x "i.) ) 



where the gain is given by 



k 








( 3 . 15 ) 




Figure 3.5. Noise Filtering by Use of Local Statistics 



Figure 3.5 shows a block diagram for the implementation of 
the algorithm. Starting with the degraded pixel Zjj,the local mean 

is evaluated and subtracted to obtain the local contrast. Then, 
both Qjj and kjj can be evaluated, and the value of k thus 

obtained is multiplied by the local contrast. This in turn is added 
to the local mean to get the restored pixel Xjj. 

One important parameter is the size of the window over 
which the local mean is estimated. If the window is too small, the 
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noise filtering algorithm is not effective. If the window is too large, 
subtle details of the image will be lost in the filtering process. 

D. MEDIAN FILTERING 

Median Filtering is a non-linear signal processing technique 
originally developed by J.Tukey [Ref. 7], that is useful for noise 
suppression in images. In one-dimensional form the median filter 
consists of a sliding window encompassing an odd number of pixels. 
The center pixel in the window is replaced by the median of the 
pixels in the window. The median of a discrete sequence a^,a 2 , 

a N for M odd is that member of the sequence for which (M-l)/2 
elements are smaller or equal in value, and (M-l)/2 elements are 
larger or equal in value [Ref. 2]. For example, if the ordered 
values of the pixels within a window are 80, 90, 200, 110, 120, 
the center pixel would be replaced by the value 110, which is the 
middle point of the sorted sequence 80,90,110,120,200. 

The concept of median filter is easily extended to two 
dimensions by utilizing a two-dimensional window of some desired 
shape such as a rectangle or a discrete approximation to a circle. 
A two-dimensional LXL median filter will provide a greater degree 
of noise suppression than sequential horizontal and vertical 
processing with LX1 median filters, but two-dimensional processing 
also results in greater signal suppression. 

The median filter has been found to be more effective than a 
linear filter for smoothing images with spiky noise degradations 
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because of extrema rejection by the median [Ref. 8] . Furthermore, 
the median filter preserves monotonic step edges, that is, it does 
not blur sharp edges as a linear low-pass filter would. 

Another interesting aspect of median filtering that has been 
studied [Ref. 9], is that of the convergence of a filtered image to 
what has been called a root signal . A root is a signal which is 
invariant under filtering by a particular median filter. The 
technique consists in making succesive passes of a noise-degraded 
image thru a median filter until an image corresponding to a root 
signal is achieved. Definition of a root signal follows, where a 
window width of 2N+1 is used. This definition uses the following 
ideas: 

1) A constant neighborhood is a region of at least N+l consecutive 
identically valued samples. 

2) An edge is an increasing or decreasing sequence of samples 
which is immediately preceded and followed by constant 
neighborhoods. An edge cannot contain any constant 
neighborhood . 

3) An impulse is a sequence of at most N consecutive samples 
whose values are different from those of the two surrounding 
regions; the two surrounding regions are identically valued 
constant neighborhoods. 

A signal is a root if and only if it contains only edges 
alternating with constant neighborhoods [Ref. 10]. The median 
filtering preserves edges and constant neighborhoods but eliminates 
impulses. The technique is based on the idea that passing a noise 
corrupted root signal through the median filter a sufficient number 
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of times will produce another root signal, which is usually “close” 
to the original root signal. In reference [Ref. 10], the concept of 
“closeness” and a measure of the number of runs to be used are 
calculated as a function of the image and window sizes. 

In this thesis, a median filter included in the Spider library 
[Ref. 11] is used, which is the implementation of the algorithm 
suggested by Huang, Yang and Tang [Ref. 12]. It consists 
primarily of an efficient way of updating the histogram each time 
that a given pixel is replaced by the median within a window. The 
algorithm is described briefly for a 3X3 window: 

Step 1: Obtain the histogram in the first window [Figure 3.6 

(a)] and find the median MDN. Next, count the pixels 
with values smaller than the median of the window. 
LTMDN is that number. 

Step 2: Shift the window right by one pixel [Figure 3.6(b)] 

and update the histogram and LTMDN. First decre- 
ment the counted values of the histogram equivalent 
to gray level values {gJ (a), gl(d), gl(g) ) of pixels a, d 
and g. That is: 

HIST [gl(a) ] = HIST [gl(a) ]-l 
HIST [gl(d) ] = msi[gl(d) ]-l 

m^[gl(g) ] = hist[^ ]-l 
LTMDN is updated as follows: 

IF gl(a) < MDN , LTMDN = LTMDN- 1 
if gl(d) < MDN , LTMDN = LTMDN-1 
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IF gl(g) < MDN , LTMDN = LTMDN-1 
In the same way for pixels and i', 

hist [gJfc'J ] = msi[gJfc'J ]+l 
msi[gJf/'J ] = H\si[gJ(f) ]+l 

msi[gJ0'J ] = msi[glO'J ]+l 

IF gl(c') < MDN , LTMDN = LTMDN+1 

IF gl(f') < MDN , LTMDN = LTMDN+1 

IF gl(i') < MDN , LTMDN = LTMDN+1 




Step 3: Update MDN in the previous window, and obtain a 

median in the new window. Let ITH= (number of 
pixels in the window) /2 
- If LTMDN > ITH, then 

LTMDN = LTMDN - HIST [MDN] 
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MDN = MDN-1 



Repeat this until LTMDN < 1TH is obtained 
- If LTMDN £ ITH, then 
MDN = MDN+1 

LTMDN = LTMDN + hist[MDN] 

Repeat this until LTMDN + hist[MDN+1] > ITH is 
obtained. This MDN is the median for the current 
window and substitute it in the output image. 

Step 4: End when one line is finished. Go to Step 2. 
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IV. EXPERIMENTAL RESULTS 



In this chapter, experimental results for different combinations 
of the algorithms described are presented and discussed. Basically, 
the adaptive filtering algorithm [Ref. 4] was used for contrast 
enhancement, as described in Chapter II. As previously stated, 
this algorithm yielded very good contrast enhancement, but also 
tended to accentuate the noise. Thus, methods for noise reduction, 
such as short space spatial subtraction [Ref. 5], noise filtering by 
use of local statistics [Ref. 6], and median filtering are used, as 
post-processors to the contrast enhancement operation. The results 
are shown in next sections. Another technique for contrast 
manipulation, contrast enhancement by use of local statistics [Ref. 
6] was also used, as an alternative to the adaptive filtering 
algorithm. These results are also given and compared to those 
previously obtained. 

A. PROCESS TYPE 1 

This section shows a combination of the adaptive filtering and 
short space spatial subtraction algorithms. In Figure 2.2 two 
degraded images are shown. The images are characterized by poor 
contrast due to the turbid water viewing conditions. These images 
were processed for contrast enhancement using the adaptive 
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filtering technique, but the noisy background was heavily 
accentuated as shown in Figure 2.4. Then the spectral subtraction 
algorithm was used to smooth the noise resulting in the images 
shown in Figure 4 . 1 . Inspection of these images shows that this 
particular algorithm was succesful in smoothing the background 
noise present in the contrast-enhanced images, and, at the same 
time, most of the information is preserved. Particularly note, in 
Figure 4.1(a), how details in the beam-shaped object in the center 
of the image are brought back, as well as those of the “fish” in 
the lower right corner of the image. The same observation can be 
made about Figure 4.1(b), in which the noise reduction was 
achieved, at the cost of some signal degradation, in this case. 

This algorithm requires the power spectral density of the noise 
to be estimated. In this case, this was accomplished by taking a 
portion of the noisy background, without any objects and using it 
to estimate the spectrum. Figure 4.2 shows a representation of the 
estimated spectrum of the noise density thus obtained. In the 
figure, the magnitude of the power spectrum is represented as the 
intensity of a 2-D image, each of the image dimensions being the 
corresponding spatial frequency coordinates. 

For the short space implementation, it is required to process 
the degraded image by dividing it into subimages and processing 
each subimage separately. In dividing the degraded image, the size 
of the subimage must be such that the image luminance in it can 
be approximated by a stationary random process. In addition, the 
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Figure 4 . 




(a) 




(b) 



Images in Figure 2.4 Processed with the Short 
Space Spatial Subtraction Algorithm 
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Figure 4.2. Power Spectrum of the Noise in Figure 2.4(a) 

window function must be a smooth function in order to avoid 
possible discontinuities that may appear at the subimage boundaries 
in the processed image. Thus, two types of smooth functions were 
used- separable 2-D triangular and Hanning windows-of size 32x32 
pixels, overlapped with its neighboring window by half the window 
duration in each dimension. It was noted that there was no 
appreciable difference in using either window -triangular or 
Hanning. 

The parameter k, appearing in Equation 3.1 is a scaling factor 
that normalizes the power spectral densities, and was evaluated as 

OO 00 

k = 2 (4.1) 

1 ,— 1 ,— 

where w(li,l 2 ) is the window function used. 
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B. PROCESS TYPE 2 



In this section, a combination of the adaptive filtering for 
Image enhancement and noise filtering by use of local statistics is 
presented. The images in Figure 2.4 were processed with the noise 
filtering algorithm, using different values for the noise variance. 
Based in the data used for generating the noise power spectrum 
shown in Figure 4.1, a noise variance of 60 was estimated, for the 
image in Figure 2.4(a). For the image in Figure 2.4(b), a noise 
variance of 45 was estimated. The algorithm was found to be quite 
sensitive to the size of the window used to evaluate the local 
statistics, the a posteriori mean and variance. Window sizes of 
8X8, 16X16, and 24X24 pixels were used without noticing any 
appreciable improvement in the processed image. Figure 4.3 shows 
the results of processing the images in Figure 2.4 using a window 
size of 32X32 and the noise variances mentioned above. Note that 
a relative reduction in the background noise was achieved, with the 
processed image resembling the effect of low-pass filtering in the 
background, but without blurring the edges of the objects in the 
images. When larger sizes for the window were used, some 
distortion in the images begun to appear. 

C. PROCESS TYPE 3 

In this section, the result of applying a median filter to a 
contrast-enhanced image, using the adaptive filtering algorithm, is 
presented. The image in Figure 2.4(a) was processed with the 
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(a) 





(b) 



Figure 4 3. Images in Figure 2.4 Processed with Noise Filtering 
by Use of Local Statistics Algorithm 
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median filter available in the Spider Library, using a window size 
of 5X5, and making succesive iterations of the processed image, 
until appreciable image degradation start to show up. The image 
obtained after 8 iterations is shown in Figure 4 . 4 . It can be 
observed that some noise smoothing was obtained. Even when the 
resulting image is not noise-free, the remaining noise is of much 
lower spatial frequency than that of the image before the 
median-filtering, and thus seems less objectionable. 

Detailed characteristics of the signal are preserved and some 
details have actually been enhanced. This is true in the case of the 
fish in the lower right corner of the image. A comparison with the 
use of 3X3 window for the filter operation led to the conclusion 
that the smaller window size gives slightly more fidelity in the 
signal preservation. However, since convergence for the 3X3 
window is much slower, this choice is undesirable. A reduction in 
the dynamic range of the median-filtered image can also be 
observed, which presents an effect similar to that of low-pass 
filtering in the background region of the image. However, this 
procedure does not have the undesirable characteristic tendency of 
the low-pass process to blur the edges of the signal. 
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Figure 4.4. Median Filtering of Image in Figure 2.4(a) 

D. PROCESS TYPE 4 

In this section the application of the contrast enhancement by 
use of local statistics algorithm is presented. In Figures 4.5 and 
4.6 four low-contrast noise-degraded images are shown. All of 
them are 512 X 512 pixels in size and each pixel is represented by 
8 bits. The images in Figures 4.5 (a),(b) and 4.6(a) are 

characterized by a relatively dark background with a diffused 
effect around the objects in the pictures resulting in low contrast. 
The image in Figure 4.6(b) has an almost uniform brightness 
throughout the picture, which makes it difficult to distinguish 
between the fish and beam shapes from the uniform background. 

Figure 4.7(a) corresponds to the image in Figure 4.5(a) 
processed with the contrast enhancement algorithm, using a value 
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of 4 for the gain k of Equation 3.10. (Recall that k is the ratio of 
the local standard deviation of the processed image to the original 
standard deviation) . Note the enhancement in contrast achieved, 
which is particularly evident in the cable connecting the two 
objects in the lower part of the image. 

The image in Figure 4.5(b) was processed using a value of 
k=6, with results shown in Figure 4.7(b). The higher value was 
used to try to enhance the fish shape in the lower right part of 
the image. This enhancement was achieved to some extent. 

Figure 4.8(a) shows the result of processing the image in 
Figure 4.6(a). In this case, a value for k of 4 was again used, 
producing an image that is a much crisper version of the 
corresponding degraded image. Also note that details, such as the 
particular shape of the object in the picture, are enhanced. The 
noise in the background is not greater than that in the original, 
which means that the desired result of contrast enhancement was 
achieved without a significant noise increase. 

The results of processing the image in Figure 4.6(b) are shown 
in Figure 4.8(b). In this case a value of k=5 was used, and the 
contrast of the two objects in the picture are increased with 
practically no increase in the noise. 
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(a) 






(b) 



Figure 4.5. Original Degraded Images 
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(b) 

Figure 4 6. Original Degraded Images 
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(b) 



Figure 4.7. Images in Figure 4 5 Processed with the Contrast 
Enhancement by Use of Local Statistics Algorithm 
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(b) 



Figure 4.8 Images in Figure 4.6 Processed with The Contrast 
Enhancement by Use of Local Statistics Algorithm 
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V. CONCLUSIONS 



Restoration of images degraded by turbid water was considered 
in this thesis. Previous work on the subject included the use of an 
adaptive filtering algorithm for enhancing the local contrast. This 
algorithm was applied to images characterized by poor contrast due 
to the turbid water viewing conditions. The technique yielded good 
contrast enhancement, but tended to accentuate the noise. 

In our work some techniques for noise smoothing were used as 
post-processors to the contrast enhancement operation. 

The specific techniques used were short space spectral 
subtraction, noise filtering by use of local statistics, and median 
filtering. The spectral subtraction algorithm yields very good noise 
reduction, but tends to introduce some signal degradation. The 
filtering technique that uses local statistics -local mean and 
variance- produced noise smoothing in the background. This 
technique was found to be very sensitive to the size of the window 
in which the local statistics were evaluated. Median Filtering did 
not eliminate the noise completely, but the remaining noise was of 
much lower spatial frequency than the original noise and thus 
seemed less objectionable. Better noise filtering was achieved when 
the degraded image was median-filtered iteratively until some 
signal degradation began to appear. 
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As an alternative to the adaptive filtering algorithm, an 
algorithm for contrast enhancement by use of local statistics was 
used. This particular technique yielded very good results without 
significant noise increase. 

All processes implemented in this thesis proved to be quite 
computationally intensive. This was especially true for the 
techniques that require the evaluation of local statistics. The 
ultimate goal for the problem addressed in this and previous work, 
is to perform the enhancement in real time. This suggests a 
natural continuation for the work done so far on this topic, in 
order to implement the algorithms for enhancement in real time. 
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APPENDIX 



COMPUTER PROGRAMS 



PR06RAM SPATSUB 



THIS PROGRAM IS THE IMPLEMENTATION OF THE SHORT SPACE 
SPECTRAL SUBTRACTION ALGORITHM 



C 

BYTE BYTEIMG(512.5I2) 

CHARACTER* 1 COEN 

CHARACTER* 40 INFILENAME.OUTFILENAME 
INTEGER* 4 IMGBUFF(5I2.512). INTARRAY(512.5I2) 



C * * * INPUT FILENAMES : INPUT & OUTPUT * * * 

TYPE 10 

10 FORMAT C INPUT FILENAME => '.$) 

ACCEPT 15. INFILENAME 
15 FORMAT (A40) 

C ACCEPT 67. IWY 

67 FORMAT (12) 

TYPE 68 

68 FORMAT (' COLUMN SIZE OF WINDOW => '.$) 
IWX =16 

C ACCEPT 69. IWX 

69 FORMAT (12) 



IERR = 0 

C * * * INPUT ORIGINAL IMAGE * * * 

CALL INPUTIMG (INFILENAME .BYTEIMG, IMGX. IMGY) 

C * * * CHANGE IMAGE FROM BYTE TYPE TO INTEGER TYPE * * * 

CALL BYTE_TOJNTEGER (BYTEIMG. IMGBUFF. IMGX. IMGY. IERR) 
IF (IERRR .EQ. 1) GO TO 100 

C * * « PERFORM THE SPATIAL SUBSTRACTION * * « 
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CALL CORTO (IMGBUFF, INTARRAY, IM6X. IHGY, IWY, IERR) 



IF (IERR .EQ. 1 ) GO TO 100 

C * * * CHANGE IMAGE FROM INTEGER TYPE TO BYTE TYPE * * * 

CALL INTEGER_TO_BYTE (INTARRAY, BYTEIMG, IMGX, IMGY.IERR) 

IF (IERR .EQ. 1) GO TO 100 

C * * * OUTPUT THE PROCESSED IMAGE * * * 

CALL OUTPUTING (OUTFILENAME, BYTEIMG, IMGX, IMGY) 

STOP 
CONTINUE 

TYPE !! FORCED TO EXIT !! ' 

STOP 

END 

SUBROUTINE CORTO (ENTRADA, SALIDA.lt16X.lt16Y.IWY .IERR) 

THIS SUBROUTINE PERFORMS THE SPATIAL SUBTRACTION ALGORITHN 
USING THE SHORT SPACE IMPLEMENTATION TECHNIQUE FOR IMAGE 
RESTORATION 



INTEGERS ENTRADA (IMGX, IMGY), SALIDA(IMGX,IMGY), INTE(32,32) 

REAL *4 AR(32,32),AI(32,32),BR(32,32),BI(32,32), 

♦ ST(64), CT(64), LBR(64), BRR(32.32), WIND0W(32,32), 

+ RENTRADA(5 1 2,5 1 2), IMGWIN(32,32), IARRAY(32,32). MAG(32,32), 

♦ PHAS(32,32),PSD(32,32) 

DO I - 1 .IMGX 
DO J - 1 .IMGY 

RENTRADA (l,J) = FLOATJ(ENTRADA(l,J)) 

END DO 
END DO 

TYPE •; INTEGER IMAGE, CONVERTED TO REAL DATA' 

C CALL MCHECK(RENTRADAJMGX) 

C * * - BUILD UP THE WINDOW FUNCTION * « » 

CALL BLTWINDOW (16, WINDOW) 

C TYPE * .WINDOW FUNCTION BUILT' 

C CALL MCHECK(WIND0W,32) 
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AK * 0.0 



001 = 1.32 
DO J *1.32 

AK * AK+(WIND0W(I.J)««2 

END DO 
END DO 



TYPE *,'AK*>... '.AK 



C * * « INITIALIZE OUTPUT ARRAY « « * 

DO L * 1 .IMG 
DO M - 1.IMGY 

SAL IDA (L.M) = 0 

END DO 
END DO 

C TYPE *. 'OUTPUT ARRAY INITIALIZED' 

C CALL MCHECK(SALIDA.IMGX) 

C *«* ESTIMATE POWER SPECTRUM DENSITY Of THE NOISE * * « 

CALL ESTINOISE(RENTRADA,PSD) 

C TYPE * .POWER SPECTRUM DENSITY Of THE NOISE ESTIMATED' 

C CALL MCHECKCPSD.32) 

C MAIN PROCESS SECTION «•* 

C INSERT LIMITS HERE 

DO I - 1,31 
DO J = 1.31 

LL * ((H)* 16H1 
MM * ((J-1)* 16H1 
DOL * 1.32 
DOM - 1,32 

IA * LL+L-1 
B * MM+M-1 

C IF(IA .GE. 16 .AND IA .LE. 112 .AND. 

C + IB .GE. 16 .AND. IB .LE. 112) THEN 

IM6WIN(L.M)-RENTRADA(IA 1 IB) 

I ARRA Y(L ,M )-WINDOW(L ,M )* IMGWIN(L ,M ) 

ELSE IF(IA GE.16 .AND. IA .LE. 1 12) THEN 

lARRAY(L.M)* 0.5*RENTRADA(IA.IB) 
ELSE IF(IB .GE. 16 .AND. IB LE. 1 12) THEN 

IARRAY(L,MH).5*RENTRADA(IA 1 IB) 
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ELSE 

I ARRAY (L,M)=RENTRADA(IA,IB) 



END IF 

END DO 
END DO 

C TYPE * .DATA ARRAY HAS BEEN WINDOWED' 

C CALL MCHECK (IARRAY.32) 

C * * * TRANSFORM DATA ARRAY INTO COMPLEX FORMAT * * * 



DO L= 1 ,32 
DO M=l ,32 

AI(L,M) = 0. 

AR(L,M) = lARRAY(L.M) 

END DO 
END DO 



C TYPE*, 'CHECK REAL PART' 

C CALL MCHECK (AR.32) 

C TYPE *, 'CHECK IMAGINARY PART (ZEROS)' 

C CALL MCHECK (AI.32) 

C * * * EVALUATE 2-D FFT OF DATA IMAGE * * * 

CALL FFTS2(AR,AI,BR,BI,32,32,ST,CT,LBR,64,2,JERR) 

C CALL MCHECK (BR.32) 

DO IX = 1,32 
DO IY = 1,32 

BRRUX.IYKABS (CMPLX(BR(IX,IY),BI(IX.IY))) 

END DO 
END DO 

C TYPE *,TFT DONE .CHECK MAGNITUDE' 

C CALL MCHECK (BRR.32) 

C CALL DISPECT (BRR) 

C * * * PERf ORFM SPATIAL SUBSTRACTION * * * 

ALFA = 1 .5 

C PSD = 40.0 

DO IX = 1,32 
DO IY = 1 .32 

BR(IX.IY) = ((BRR(IX,IY))**2) - ((ALFA/AK) * ((PSD(IX.IY) **2))) 
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ELSE 

MAG(IX.IY) = SQRT(BRdX.IY)) 

END IF 

END DO 
END DO 

TYPE -.SPATIAL SUBT.DONE CHECK RESULT' 

CALL MCHECK (MAG) 

»»• EVAL. MAGNITUDE OF INV, FFT--» 

DO IX= 1 ,32 
DO IY= 1,32 

BR(IX.IY) = MAG(IX.IY) 1 * COS(PHAS(IX,IY) 

Bl (IX, IY) = MAG(IX.IY)* SIN(PHAS(IX,I Y) 

END DO 
END DO 

CALL FFTS2(BR.BI .AR.AI ,32.32,ST.CT,LBR,64,-2,JERR) 

TYPE", ‘INV. FFT DONE. CHECK REAL PART' 

CALL MCHECK(AR) 

TYPE -.'CHECK IMAGINARY PART (ZEROS) 1 
CALL MCHECK(AI) 

- » » BACK TO INTEGER FORMAT * - - 

DO IX = 1,32 
DO IY = 1,32 

INTE(IX.IY) = JINT(AR(IX.IY)) 

END DO 
END DO 

TYPE- .'INTEGER CONVERSION DONE, CHECK IT' 

CALL MCHECK (INTE) 

- - » BUILD UP THE OUTPUT ARRAY - - - 

DO IX = 1 ,32 
DO IY = 1.32 

LL =((M)-16H 
MM((J-I )- 16}+ 1 
IXX =LL+IX-1 
IYY =MM+IY-1 

SALIDA (IXX, IYY) = SALID AdXX.I YYHINTEdX.I Y) 
TYPE-, 'IXX AND IYY =>',IXX,IYY 



oooooo oo ooooo 



END DO 
END DO 



END DO !! END DO J 
END DO !! END DO I 

100 TYPE*,T1AIN PROCESS DONE' 

C CALL TRUNCATE(SALIDA,IH6S,IM6Y) 

CALL SCALE_ADJUST(SALIDA.IMGX,IMGY, 0.0, 255.0) 
TYPE * .'SPATIAL PROCESS DONE' 

RETURN 

END 



SUBROUTINE BLTWINDOW(IWY.NIN2) 

REAL "4 NK32). NIN2C32.32), N2(32,32) 

L =0 

PI = 4*ATAN (1.) 

DO I + 1 .32 

II =1-1 

N 1 (I ) = 0.5 - (0.5** COS(2" PI *• 11/3 1 )) 
N2(l) * N1(l) 

TYPE *, N 1(1), N2(l) 

IF (I ,GT. 16) THEN 
L = L+1 
N 1 (I ) = I - (2*L) 

ELSE 

N 1 (I) = I 
END IF 

END DO 
L =0 

DO I = I . 32 
DO J = 1,32 

IF (J .GT. 16) THEN 
L =L +1 

NIN2(U) = Nl(l) * (J-2* 1 ) 

ELSE 

NIN20.J) = Nl(l) * J 

END IF 

NIN2(U) = N1(I)-N2(J) 

C TYPE *, NIN20.J) 

END DO 
L -0 
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END DO 



RETURN 

END 



SUBROUTINE ESTINOISE (RENTRADA.PSD) 



THIS ROUTINE ESTIMATES THE POWER SPECTRAL DENSITY OF 
THE NOISE IN AN IMAGE. BY TAKING THE MAGNITUDE (SQUARED) 

Of THE DISCRETE SPACE FOURIER TRANSFORM Of A (NOISY) PART 
OF THE BACKGROUND. 

C TO CHOOSE APPROPIATE REGION. SET PARAMETERS I AND J 

C ACCORDINGLY. 

REALM RENTRADA(5 1 2.5 12), PSD(32.32). AR(32.32). AK32.32). 
BR(32,32), BK32.32). ST(32.32), CT(32.32), LBR(64) 

C "*• TAKE A PORTION Of THE IMA6E BACKGROUND « « « 



J-80 

K-80 

DO L-1 .32 
DO M= 1 .32 

LL-L+J 

MM=M+K 

AR(L.M)-RENTRADA(LL.MM) 

END DO 
END DO 

C • » • INITIALIZE IMAGINARY PART Of DATA ■ • * 

DO L-1. 32 
DO M-1,32 

AKL.MH). 

END DO 
END DO 

C ■ • « EVALUATE MAGNI TUDE OF 2-0 FFT » - * 

CALL FFTS2(AR.AI.BR.BI.32.32.ST,CT.LBR,64,2.JERR) 

DO L-1. 32 
DO M-1,32 

PSD(L ,MK ABS(CMPLX(BR(L .M ).BI (L ,M ))) 

END DO 
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END DO 



RETURN 

END 



SUBROUTINE INPUTIf16(INFILENAhE.BYTEItt6.IM6X.IM6Y) 

C READ IN INPUT IMAGE DATA IN BYTE TYPE 

C ASSUME INPUT DATA FILE IS A DIRECT ACCESS FILE 

BYTE BYTEIMG(IM6Y,IMGX) 

CHARACTER* 40 INFILENAME 



0PEN(UNIT=2, FILE=INFILENAME, STATUS='OLD' ( ACCESS-DIRECT* . 
+ RECORDTYPE=TIXED') 



DO M.IMGY 

READ(2'I) (BYTEIM6(I, J). J=1, IMGX) 

END DO 

TYPE * , ‘ * INPUT IMAGE HAS BEEN READ ' 

CL0SE(2) 

RETURN 

END 



SUBROUTINE BYTE_T0_INTE6ER (BYTEDATA. INTDATA. IX. IY. IERR) 

THIS PROGRAM CHANGE DATA (USUALLY 2D IMAGE DATA) IN THE BYTE DATA 
TYPE INTO THE INTEGER TYPE 



BYTE BYTEDATA (IY. IX) 
INTEGER-4 INTDATA (IY. IX) 



DO I = I, IY 
DO J = I. IX 

IF (BYTEDATA (I.J) .GE. -128 
+ .AND. BYTEDATA (U) IT. 0) THEN 
INTDATAO.J) = BYTEDATA (I.J) + 256 
ELSE IF (BYTEDATA(U) .GE.O 
+ .AND. BYTEDATA(U) ,LE. 127) THEN 
INTDATA(U) = BYTEDATA(U) 

ELSE 
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TYPE* ,'TJDATA OUT OF RANGE! R CM, COL. DATA => 

+ I. J. BYTEDATA(I.J) 

IERR = 1 
GOTO 100 
ENDIF 

END DO 
END DO 

TYPE*. ' BYTE TO INTEGER CONVERSION HAS BEEN DONE - 

100 CONTINUE 

RETURN 

END 

SUBROUTINE INTE6ER_T0JBYTE(INTDATA, BYTEDATA. IX. IY. IERR) 

THIS PROGRAM CHANGE DATA (USUALLY 2D IMAGE DATA) IN THE BYTE DATA 
TYPE INTO THE INTEGER TYPE 



BYTE BYTEDATAOX. IY) 
INTEGER* A INTDATAOX, IY) 



DO 1=1 .IX 
DO J=1 ,IY 



IF (INTDATA(U) .LT. 0 .OR. INTDATA(U) .GT. 255) THEN TYPE*. - S 
DATA OUTOF RANGE! ROW.COL.DATA => 

+ I. J, INTDATA(I.J) 

IERR = 1 
GOTO 100 

ELSE IF (INTDATAO.J) .IE. 127) THEN 
BYTEDATA(I.J) * INTDATAU.J) 

ELSE 

BYTEDATA(U) = INTDATA(U) - 256 
ENDIF 

END DO 
END DO 

TYPE *.' INTEGER TO BYTE CONVERSION HAS BEEN DONE - 

100 CONTINUE 
RETURN 
END 
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PROGRAM LOC-STAT 



C THIS PROGRAM IS THE IMPLEMENTATION Of THE IMAGE ENHANCEMENT AND 

C NOISE FILTERING ALGORITHM 

C 
C 

BYTE BYTEIMG (512,512) 

CHARACTER* 1 COEN 

CHARACTER* 40 INFILENAME, OUTFILENAME 

INTEGER* 2 IMGBUFF(5 1 2,5 1 2), LAVEIMG(5 1 2.5 1 2), L0CMEAN(512,512), 

- INTARRAYC5 1 2,5 1 2) 



C * * * INPUT FILENAMES : INPUT & OUTPUT * * * 

TYPE 10 

10 FORMAT C INPUT FILENAME -> ',$) 

ACCEPT 15, INFILENAME 

15 FORMAT (A40) 

TYPE 20 

20 FORMAT (' OUTPUT FILENAME -> \X) 

ACCEPT 25, OUTFILENAME 

25 FORMAT (A40) 

C * * * INPUT ORIGINAL IMAGE SIZE * * * 

TYPE 50 

50 FORMAT C ROW SIZE Of IMAGE => ',$) 

ACCEPT 55. IMGY 

55 FORMAT (14) 

TYPE 60 

60 FORMAT C COLUMN SIZE Of IMAGE => ',$) 

ACCEPT 65, IMGX 

65 FORMAT 914) 

C *** INPUT WINDOW SIZE *** 

TYPE 66 

66 FORMAT CROW SIZE Of WINDOW ->',$) 

ACCEPT 67, IWY 

67 FORMAT (12) 

TYPE 68 

68 FORMAT (’ COLUMN SIZE Of WINDOW => ',$) 

ACCEPT 69, IWX 

69 FORMAT (12) 

C * * - SELECT CONTRAST ENHANCEMENT OR NOISE FILTERING * * * 

TYPE 70 

70 FORMAT ( 'CONTRAST/NOISEFIL ? (C/N) »> \ $) 
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ACCEPT 75. CE 
75 FORMAT (A 1) 



IERR = 0 

C *** INPUT ORIGINAL IMAGE * * - 

CALL INPUTIM6 (INFILENAME. BYTEIMG, IMGX. IM6Y) 

C * * * CHANGE IMAGE FROM BYTE TYPE INTEGER TYPE * * * 

CALL BYTE-TOJNTEGER ( BYTEIMG. IMGBUFF. IMGX. IMGY, IERR) 

IF (IERR EQ.DGOTO 100 

C * * * CALCULATE LOCAL MEAN ARRAY * * * 

CALL LOCAL_MEAN(IMGBUFF, LOCMEAN, IMGX. IMGY. IWX. IWY, IERR) 
IFOERR EQ.1 ) GOTO 100 

IF (CE EQ. 'C .OR. CE .EQ. c ) THEN 

C - * * INPUT K = STANDARD DEVIATION RATIO « - * 

TYPE 80 

80 FORMAT (' STAND. DEV. RATIO K->\ $) 

ACCEPT 85. K 
85 FORMAT (M) 

CALL ENHACONT (IMGBUFF. LOCMEAN. INTARRAY. IMGX, IMGY. K. IERR) 
ELSE 

C " * * INPUT NV = NOISE VARIANCE * * * 

TYPE 90 

90 FORMAT C NOISE VARIANCE NV =>',$) 

ACCEPT 95 

95 FORMAT (14) 

CALL NOISE F ILT ( IMGBUFF. LOCMEAN. INTARRAY. NV. IMGS. IMGY. IERR) 

ENDIF 

SCLMIN - 0.0 
SCLMAX = 255.0 

C CALL SCALE_ADJUST (INTARRAY, IM6X, IMGY, SCLMAX. SCLMIN) 

CALL TRUNCATE (INTARRAY, IMGX, IMGY) 

C * » * CHANGE IMAGE FROM INTEGER TYPE TO BYTE TYPE * « * 

CALL INTEGER_TO_BYTE ( INTARRAY, BYTEIMG. IMGX. IMGY. IERR ) 

IF (IERR .EQ. DGOTO 100 

C « « * OUTPUT THE PROCESSED IMAGE « * - 
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CALL OUTPUTIMG ( OUTFILENAME, BYTEIMG, IMGX, IMGY ) 



STOP 

100 CONTINUE 

TYPE *.' !! FORCED TO EXIT !! ‘ 
STOP 

END 



SUBROUTINE ENHACONT (IM6BUFF.L0CMEAN.INTARRAY.in6X,IM6Y,K.IERR) 



THIS SUBROUTINE IMPLEMENTS THE CONTRAST ENHANCEMENT BY USE 
OF LOCAL STATISTICS ALGORITHM. 



C 

INTEGER* 2 IMGBUFF(IMGYJMGX) . LOCMEAN(IMGYJMGX) , 
INTARRAY(IMGY,IMGX) 

C * * * EVALUATE LOCAL CONTRAST AND MODIFY IT WITH K * * * 

DO I - 1 IMGY 
DO J » I, IMGX 

INTARRAYCU) » K*(IMGBUFF(I,J) - LOCMEAN(U)) 

END DO 
END DO 

C * * * ADD MODIFIED LOC . CONTRAST AND LOCAL MEAN * * * 

DO I * UMGY 
DO J = I .IMGX 

INTARRAY (U) - INTARRAY (l,J) + LOCMEAN (I.J) 

END DO 
END DO 

C * * * CORRECT ANY OVERFLOW * * » 

DO I - UMGY 
DO J = IJM6X 

INTARRAY(U) = K*(IM6BUFF(I.J) - LOCMEAN(U) ) 

END DO 
END DO 

C * » » ADD MODIFIED LOC. CONTRAST AND LOCAL MEAN * * * 

DO I = UMGY 
DO J = I.IMGX 

INTARRAY (I.J) = INTARRAY (l,J) + LOCMEAN (I.J) 

END DO 
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END DO 



»*« CORRECT ANY OVERFLOW * * * 

SCLMIN = 0.0 
SCLMAX = 255.0 

CALL SCALE_ADJUST (INTARRAY. IMGX. IMGY. SCLMAX, SCLMIN) 

RETURN 

END 



SUBROUTINE N0ISE_FILT(lf16BUFF.L0CMEAN.INTARRAY,NV.IM6X, 
+ IM6Y.IERR) 



THIS SUBROUTINE IS THE IMPLEMENTATION OF THE NOISE FILTERING 
ALGORITHM, CALLED BY THE NOISE FILTERING AND CONTRAST ENHANCEMENT 
BY USE OF LOCAL STATISTICS PROGRAM. 

REFERENCES: 



INTEGER* 2 IMGBUFF(IMGY.IMGX), LOCMEANOMGY.IMGX), 
INTARRAY(IMGY.IMGX) 

REAL * 8 AK( 128, 128) .0(128.128) 

* * » EVALUATE THE VARIANCE OF THE IMPUT IMAGE * * * 

TYPE * .‘CHECK IMGBUFF ARRAY' 

CALL MCHECK(IMGBUFF) 



C 

16 

C 

17 



TYPE * , ‘CHECK LOCMEAN ARRAY ‘ 
CALL MCHECK(LOCMEAN) 

DO I = 1 .IMGY 
DO J = 1 ,IMGX 

IDIF = IMGBUFF(I.J) - LOCMEANO.J) 

TYPE 16, IDIF 

FORMAT C DIF => 110) 

(XU) = (IDIF* * 2) - NV 
TYPE 17, (XU) 

FORMATC Q => ‘.F8.2) 

END DO 
END DO 



C 

C 



TYPE *,' CHECK THE Q ARRAY ' 
CALL MRCHECK(Q) 

* * « EVALUATE THE GAIN AK * * « 



DO 1= 1 .IMGY 
DO J = 1 , IMGX 
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DEN = 0(1. J) + NV 
IF ( DEN ,EQ. 0.0) THEN 
DEN = 1 .0 

ENDIF 

AKCI.J) = (0(1, J) / DEN) 

END DO 
END DO 

C *** ESTIMATE THE RESTORED IMAGE * * * 

DO I = 1. IMGY 
DO J = I , IMGX 

DIF = AK(I,J)*(IM6BUFF(I,J) -LOCMEANO.J)) 
INTARRAYO.J) = INT(LOCMEANO.J) + DIF) 

END DO 
END DO 

RETURN 

END 



SUBROUTINE SCALE-ADJUST (IM6BUFF, IM6X. IM6Y. SCLMAX. SCLMIN) 

INTEGERS IMGBUFF (IMGX, IMGY) 

MAXVAL = IMGBUFF (1,1) 

MINVAL * IMGBUFF (1,1) 

DO I =1, IMGY 
DO J = 1 , IMGX 

IF (IMGBUFF(I,J) .GT. MAXVAL) THEN 
MAXVAL = IMGBUFF(U) 



ENDIF 

IF (IMGBUFF(IJ) .IT. MINVAL) THEN 
MINVAL = IMGBUFF(U) 

ENDIF 



END DO 
END DO 

C TYPE 10. MAXVAL. MINVAL 

CIO FORMATCMAX &MIN VALUES =>' 218) 

TYPE \T1AXVAL= '.MAXVAL, T1INVAAL= '.MINVAL 

INTVAL = MAXVAL - MINVAL 

SLOPE = (SCLMAX - SCLMINVREAL (INTVAL) 
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DO I = I , IMGY 
DO J « 1, IMGX 



IF (IMGBUFF(U) .EQ. MINVAL) THEN 
IMGBUFF (l,J) - SCLMIN 

ELSE IF (IMGBUFF(U) .EQ. MAXVAL) THEN 
IMGBUFF (l,J) - SCLMAX 

ELSE 

IDIF = IMGBUFF (l,J) - MINVAL 
IMGBUFF(I.J) = INT(S) OPE* IDIF) + INTCSCLMIN 

ENDIF 

END DO 
END DO 

RETURN 

END 



SUBROUTINE TRUNCATE0M6BUFF, IM6X. IM6Y) 

C TRUNCATE IMAGE GRAY VALUES INTO 0-255 

INTEGER* 4 IMGBUFFOMGY, IMGX) 

DO I - I , IMGY 
DO J - 1, IMGX 

IF (IMGBUFF(U) IT. 0) IMGBUFF(U) = 0 
IF (IMGBUFFO.J) .GT. 255) IMGBUFF(I.J) = 255 

END DO 
END DO 

RETURN 

END 
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PR06RAM HEDIFILT 



C THIS PROGRAM IS THE IMPLEMENTATION Of THE MEDIAN 

C FILTERING IN NOISE REDUCTION 



BYTE BYTEIMG(512, 512) 

CHARACTER *40 INFILENAME, OUTFILENAME 

INTEGER* 4 IMGBUFF(512,5I2). IHST(256), INTARRAY(5I2, 512) 



C * * * INPUT FILENAMES : INPUT & OUTPUT * * « 

TYPE 10 

10 FORMAT ('INPUT FILENAME => '.$) 

ACCEPT 15, INFILENAME 
FORMAT (A40) 

TYPE 20 

20 FORMAT ('OUTPUT FILENAME => ',$) 

ACCEPT 25, OUTFILENAME 
25 FORMAT (A40) 

C * * * INPUT ORIGINAL IMAGE SIZE * * * 

C TYPE 50 

50 FORMAT CROW SIZE OF IMAGE =>',$) 

IMGY=5 1 2 

C ACCEPT 55, IMGY 

55 FORMAT (14) 

C TYPE 60 

60 FORMAT C COLUMN SIZE OF IMAGE =>',$) 

IMGX = 512 

C ACCEPT 65. IMGX 

65 FORMAT (14) 

C TYPE *, IMGX, IMGY 

IERR = 0 

C * * * INPUT ORIGINAL IMAGE “ * * 

CALL IMPUTIMG (INFILENAME .BYTEIMG, IMGX, IMGY) 

C * * * CHANGE IMAGE FROM BYTE TYPE TO INTEGER TYPE * * * 

CALL BYTE-TOJNTEGER (BYTEIMG. IMGBUFF. IMGX, IMGY, IERR) 
IF (IERR ,EQ. 1 ) GO TO 100 

C * * * PERFORM THE MEDIAN FILTERING * * * 

DO I = 1,20 

CALL MEDI (IMGBUFF, INTARRAY, IMGX, IMGY. 5, 5, IHST, 256) 
DO IX = 1 , IMGX 
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DO IY = 1, IMGY 

IMGBUFF(IX.IY) = INTARRAYOX, IY) 

END DO 
END DO 

END DO 

C * * * CHANGE IMAGE FROM INTEGER TYPE TO BYTE TYPE * » * 

C TYPE MMGX, IMGY 

CALL INTEGER_TO_BYTE (INTARRAY, BYTEIMG, IMGX, IMGY. IERR) 

IF (IERR .EQ. 1) GO TO 100 

C «** OUTPUT THE PROCESSED IMAGE * * * 

CALL OUTPUTIMG (OUTFILENAME. BYTEIMG, IMGX. IMGY) 

STOP 

100 CONTINUE 

TYPE !! FORCED TO EXIT ! !’ 

STOP 

END 



62 



LIST OF REFERENCES 



1 Gonzalez, R., Digital Image Processing . Addison-Wesley, 1977. 

2. Pratt, W K., Digital Image Processing . John Wiley & Sons, 
Inc., 1978. 

3. Chellappa, R. and Sawchuk, A. A.. Digital Image Processing 
and Analysis : Volume 1: Digital Image Processing , IEEE Compu- 
ter Society Press, 1985. 

4. Peli, T. and Lim, J. S., “Adaptive Filtering for Image En- 
hancement,” Proceedings of IEEE International Conference on 
Acoustics Speech and Signal Processing , March 30, 31 8, April 
1, 1981. 

5. Lincoln Laboratory, M. I. T. ,TN 1979-17. Image Restoration bv 
Short Space Spectral Subtraction . bv J.S. Lim, 1979. 

6. Lee,J. “Digital Image Enhancement and Noise Filtering by Use 
of Local Statistics, ” IEEE Transactions on Pattern Analysis and 
Machine Intelligence . Vol. PAMI-2, Number 2, pp 165-168,1980. 

7. Tukey,J. W. . Exploratory Data Analysis . Addison-Wesley. 1971 

8. Narendra, P. M.,“A Separable Median Filter for Image Noise 
Smoothing, ” IEEE Transactions on Pattern Analysis and 
Machine Intelligence . Vol. PAMI-3, Number 1, pp 20-29, 
January 1981. 

9. Wendt, P. D. , Coyle, E.J,,and Gallagher, N.C., “Some Conver- 
gence Properties of Median Filters, ” IEEE Transactions on 
Circuits and Systems . Vol. CAS-33, Number 3, March 1986. 

10. Gallagher, N.C. and Wise, G.L.,“A theoretical Analysis of the 
Properties of Median Filters, ” IEEE Transactions on Acoustics- 
Speech and Signal Processing , Vol. ASSP-29, pp 1136-1141, 
1981. 



63 



11. Huang, T.S., Yang, G.J., and Tang, G.Y., "A Fast Two- 

Dimensional Median Filtering Algorithm, ” IEEE Transactions on 
Acoustics. Speech. and Signal Processing . Vol. ASSP-27, Number 
1, 1979. 



64 



INITIAL DISTRIBUTION LIST 



No. Copies 

1. Defense Technical Information Center 2 

Cameron Station 

Alexandria, Virginia 22304-6145 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93943-5002 

3. Department Chairman, Code 62 1 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943 

4. Professor C. W. Therrien, Code 62Ti 3 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943 

5. Professor P. H. Moose, Code 62Me 1 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943 

6. Sr. Almirante Comandante de la Armada 2 

Comando Armada Nacional 

CAN, Avenida ElDorado 
Bogota, D.E. 

Colombia 

7 Sr. Contralmirante Director Escuela Naval 2 

Escuela Naval “Almirante Padilla” 

Cartagena 

Colombia 



65 



8 . 



1 



Mr. Richard Evans, Code 7021 
NUWES Code 7021 

Naval Undersea Weapons Engineering Station 
Keyport, Washington 98345 

9. Ms. Tami Peli 1 

Gr. 21 

M.l.T. Lincoln Laboratory 
Lexington, Massachusetts 02173-0073 

10. Professor Jae Lim . 1 

Department Elec. Eng. & Computer Science 

Room 36-653 
M.l.T. 

Cambridge, Massachusetts 02139 

11. Professor 0. B. Wilson, Code 61W1 1 

Department of Physics 

Naval Postgraduate School 
Monterey, California 93943 

12. Ms. Atsuko Keiyu , Code 62 1 

Department of Electrical and Computer Engineering 
Naval Postgraduate School 

Monterey, California 93943 

13. Teniente de Navio Jorge A. Franco L. 2 

Escuela Naval “Almirante Padilla” 

Cartagena 

Colombia 



66 



DUDLEY KNOX LIBRARY 

^ ^ SCHOOL 

I. UrLRL\ - OAii FOF-NlA 93943-BOOS 



22100G 

Thesis 

F754 Franco 

c.l Enhancement of video 

images degraded by turbid 
water. 



2?:ooc 

Thesis 

F754 Franco 

C «1 Enhancement of video 

images degraded by turbid 
water. 



